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Bi-stability in turbulent, rotating spherical Couette flow 
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Flow between concentric spheres of radius ratio 77 = ri/ro — 0.35 is studied in a 3 ni outer diameter experiment. 
We have measured the torques required to maintain constant boundary speeds as well as localized wall 
shear stress, velocity, and pressure. At low Eknian number E — 2.1 x 10~^ and modest Rossby number 
0.07 < Ro < 3.4, the resulting flow is highly turbulent, with a Reynolds number {Re = Ro/E) exceeding 
fifteen million. Several turbulent flow regimes are evident as Ro is varied for fixed E. We focus our attention 
on one flow transition in particular, between Ro =1.8 and Ro = 2.6, where the flow shows bistable behavior. 
For Ro within this range, the flow undergoes intermittent transitions between the states observed alone 
at adjacent Ro outside the switching range. The two states are clearly distinguished in all measured flow 
quantities, including a striking reduction in torque demanded from the inner sphere by the state lying at 
higher Ro. The reduced angular momentum transport appears to be associated with the development of a 
fast zonal circulation near the experiment core. The lower torque state exhibits waves, one of which is similar 
to an inertial mode known for a full sphere, and another which appears to be a strongly advected Rossby-type 
wave. These results represent a new laboratory example of the overlapping existence of distinct flow states in 
high Reynolds number flow. Turbulent multiple stability and the resilience of transport barriers associated 
with zonal flows are important topics in geophysical and astrophysical contexts. 
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I. INTRODUCTION 

An understanding of rapidly rotating turbulent flow is 
key to many problems of geophysical and astrophysical 
fluid dynamics. The large scale fluid motions of stars 
as well as planetary atmospheres, oceans, and cores are 
strongly influenced by the Coriolis accelerations arising 
from the rapid overall rotation in those systems. These 
flows can be described by the Navier-Stokes equation 
written for a frame rotating with constant angular ve- 
locity rj = flz, which in dimensionless form is 



expresses the strength of the nonlinear term relative to 
the Coriolis term 2z x u, and the Ekman number, 



E^ 



{VlP) ' 



(3) 



characterizes the importance of viscous drag relative to 
Coriolis acceleration. Throughout this paper we use, ^ = 
To — r-i, Vt — Q,o, and U = (ilj — fio)^, so that 



Ro = 






and 



du 

Ik 



i?o u • Vu + 2z X u = -Vp + E V^u. (1) 



E = 



nj^' 



(4) 



(5) 



Using U as the characteristic velocity scale, £ as the 
characteristic length scale, and the kinematic viscosity 
v, we dcflnc the relevant dimensionless parameters. The 
Rossby number, 



Ro 



U 



(2) 
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The Taylor-Proudman constraint on rotating flows of- 
ten holds approximately: overall rotation about the z- 
axis leads to a tendency toward z-independence of the 
flow for large scales or slow motions. This is useful for 
nearly steady flow where (u- V)u and du/dt in Eq. 1 are 
small with respect to 2z x u. Due to the relatively large 
size of this apparatus and practical limitations, all of our 
flows are turbulent. Thus the Taylor-Proudman theorem 
is unlikely to hold, though it is expected that significant 
anisotropy will remain. When accelerations are signifi- 
cant, retaining time dependence but neglecting the non- 
linear term in Eq. 1 yields dynamics where flow distur- 
bances can propagate via linear Coriolis-restored inertial 
waves^. In containers, wave modes arise. Modes for the 



interior of a sphere are treated by Greenspan^. Zhang 
et al.^ provide complete exphcit analytical solutions for 
the full sphere. The modes in a spherical annulus are not 
known analytically but have been investigated numeri- 
cally'^"^. Inertial waves and related Rossby waves play 
a role in the inverse energy cascade to large scales ob- 
served in rotating turbulence^. Experiments have shown 
that inertial modes are important in turbulent spherical 
Couette flow^"^°, and in more general turbulent flows 
like rotating grid turbulence^^. We expect that iner- 
tial modes have substantial influence on the character 
of the flow in rapidly rotating bounded systems for mo- 
tions with frequencies less than twice the rotation rate. 
Above that frequency, no inertial waves or modes exist. 
The question of completeness of inertial modes is still an 
open one^^, so it may not be possible to express arbitrary 
motions with frequencies below a; = 2 f^o in terms of in- 
ertial mode Fourier components. However, the nonlinear 
interactions of a sea of modes with uj < 2Qo may have a 
strong signature on a wide range of spatial and temporal 
scales. 

Turbulent flow in a rapidly rotating spherical annu- 
lus with a radius ratio 77 = n/^^o — 0.35 has poten- 
tial geophysical relevance due to geometrical similarity 
to Earth's liquid outer core. Differential rotation im- 
posed by the boundaries is at flrst glance considerably 
different from the convection, precession, and tidal forc- 
ing that may drive flows in planetary systems. However, 
convection in rotating systems tends to set up differential 
rotation^'^"^^ , as does the nonlinear interaction of inertial 
mode shear layers in spherical shells driven by precession 
or tidal deformation^^'^°. 

The study of spherical Couette flow with rapid outer 
sphere rotation has been fairly limited. Experiments and 
theoretical studies have focused more on the case with 
the outer sphere stationary^^"^''. In experiments, this 
is possibly due to of the difficulty of conducting mea- 
surements in the rotating frame. Furthermore, there is 
a complicated dependence of the observed laminar flow 
on Reynolds number and gap width even when the outer 
sphere is fixed^^'^^. Nevertheless, some work in spheri- 
cal Couette with overall rotation has been carried out. 
Hollerbach, Egbers, Futterer, and More'^° teamed ex- 
perimental observations with numerical simulations to 
study Stewartson layer instabilities in the case of small 
to modest differential rotation. They found good agree- 
ment between simulation and observed visual patterns 
in experiments for parameters for E > 8 x 10^** and 
< Ro < 0.6. Experiments by Gertsenshtein, Zhilenko, 
and Krivonosova'^^ investigated the transition to turbu- 
lence at some values of Ro, especially when Ro < 0. 
Schaeffer and Cardin'^^ studied the onset of Stewartson 
layer instabilities at E similar to our experiment at very 
low Ro where the first instability happens. 

The strongly turbulent behavior of spherical Couette 
flow has been studied mostly in hydromagnetic apparatus 
with liquid sodium as a working fluid. Sisan et al.^^ stud- 
ied instabilities of a turbulent flow of sodium that arise 



with a sufficiently strong axial magnetic field, with some 
hydrodynamic measurements to characterize the initial 
unmagnetized flow. Kelley et al7'^ inferred the flow in 
a 60 cm rapidly rotating sodium apparatus with a mag- 
netic field sensor array and a dynamically passive applied 
field. These experiments demonstrated over-refiectional 
excitation of inertial modes at comparable Ekman num- 
ber to ours and —2 < Ro < 0. Rapidly rotating spherical 
Couette flow of sodium strongly magnetized by a dipole 
permanent magnet inner sphere has been studied in the 
DTS experiment^"'^'^''^^ in Grenoble. Measurements of 
velocity, magnetic fleld, and electric potential measure- 
ments exposed a number of interesting hydromagnetic 
states in the turbulent regime^"''^'^''^"'. 

At present, numerical work on spherical Couette flow 
has not investigated the portion of the Ro, E parameter 
plane where we find turbulent bi-stability. However, nu- 
merical and theoretical work has demonstrated interest- 
ing phenomena with outer sphere rotation. Stewartson'''' , 
studied the case at infinitesimal Rossby number, deriv- 
ing the form of the free cylindrical shear layer tangent 
to the inner sphere equator that still bears his name. A 
systematic three-dimensional numerical study of instabil- 
ities of the Stewartson was undertaken by Hollerbach''®. 
Schaeffer and Cardin^-''^ used depth-averaged equations 
of motion coupled to realistic Ekman pumping to con- 
duct quasi-geostrophic simulations at low Ekman number 
for a split-sphere geometry similar to spherical Couette 
fiow. Schaeffer and Cardin''^ showed good agreement be- 
tween experiment and simulation, predicting Stewartson 
layer instabilities correctly for very low Rossby number. 
Quasi-geostrophic simulations'' found Rossby wave tur- 
bulence at higher Ro. 

Our experiments can not generally achieve Ro < 0.05 
at experimentally accessible Ekman number because of 
limitations on motor minimum speeds. The lowest possi- 
ble Rossby number increases if we rotate the system more 
slowly to raise E. It is impractical to reduce the Rossby 
number or increase the Ekman number enough to match 
the parameters of known numerical simulations. 

Guervilly and Cardin^^ performed simulations of 
spherical Couette flow with outer rotation in an inves- 
tigation of magnetic dynamo action. They performed 
fully three dimensional numerical simulations for Ekman 
number higher than E = 10^^ while achieving Rossby 
number matching some presented here. The deflnition 
of the Rossby number used by Guervilly and Cardin'''^ 
is different, -Roqc — ?/Af7/i7o ~ rjRo. They report an 
771 = 2 Rossby type wave for Ro ^ 2.9 (-Rogc = l)j 
above our first bistable range. The existence and az- 
imuthal wavenumber of this wave is partially consistent 
with our observations at that Rossby number. However, 
we observe strong turbulence, additional waves, and the 
bi-stable behavior that is the focus of this paper. These 
have not been previously reported in spherical Couette 
fiow. 

Turbulent flow transitions and multiple stability in 
very high Reynolds number flows have been reported 



in a number of systems. Several examples exist in geo- 
physics and astrophysics. The dynamo generated mag- 
netic fields of the Earth and Sun reverse polarity. Ocean 
currents, namely the Kuroshio current in the North Pa- 
cific near Japan and the Gulf Stream, both exhibit bi- 
stability in meander patterns. ^^ Polar vortices in Earth's 
stratosphere are bounded by resilient transport barriers 
much of the time, mixing to higher latitudes only inter- 
mittently. ^^"^^ This process is important to ozone deple- 
tion in polar regions and may have some similarities to 
the dynamics we observe. 

A number of laboratory flows are known to exhibit 
multi-stability and hysteresis. The mean circulation in 
turbulent thermal convection cells has been observed to 
switch direction abruptly.**^ Hysteresis in the large scale 
flow states has been seen in surface waves excited by 
turbulent swirling flows in a Taylor Couette geometry 
with a free surface.*^ Von Karman flow in a cylinder be- 
tween two independently rotating impellers has exhibited 
multi-stability and hysteresis of the mean flow despite 
extremely high fluctuation levels. ^'*"'*^ Magnetohydrody- 
namic experiments in the von Karman geometry have 
succeeded in producing dynamos that show reversals of 
the generated magnetic field. '*^'''* The L-H transition in 
turbulent tokamak plasma confinement devices involves 
the formation of a wave-driven zonal flow transport bar- 
rier that greatly enhances conflnement of the plasma. ''^'^'^ 
However, this barrier eventually breaks down in a burst 
that can damage the confinement device."''^ Understand- 
ing this particular form of turbulent multiple stability 
and its control is important issue in sustained confine- 
ment of fusion plasmas. 

Spherical Couette flow is a dynamically rich system 
in both laminar and turbulent regimes. Our appara- 
tus, initially designed, constructed, and eventually des- 
tined for magnetohydrodynamic experiments in molten 
sodium metal, presents a unique opportunity to measure 
the properties of hydrodynamic turbulence in this geom- 
etry, including transition phenomena between multiple 
turbulent states. We operate in a novel region of pa- 
rameter space, simultaneously achieving moderately high 
Rossby number and low Ekman number. This regime can 
currently only be directly accessed by experiments and 
naturally occurring flows, and has not been the focus of 
previous purely hydrodynamic studies. Furthermore, we 
are able to make quantitative measurements in the rotat- 
ing frame, something that can be quite difficult in smaller 
apparatus. 



II. APPARATUS 

The three meter apparatus allows independent rota- 
tion of the inner and outer spheres. Instrumentation 
in the rotating frame allows measurements of velocity, 
wall shear and pressure, as well as the torques required 
to maintain the boundary speeds. Fig. 1 is a schematic 
sketch of the apparatus. The stainless steel outer vessel 




FIG. 1. A schematic of the apparatus showing the inner and 
outer sphere and locations of measurement ports in the vessel 
top lid at 60 cm cyhndrical radius (1.18 n). Data from sensors 
in the ports are acquired by instrumentation, including an 
acquisition computer, bolted to the rotating lid and wirelessly 
transferred to the lab frame. Also shown is the wireless torque 
sensor on the inner shaft. 



has an inner diameter of 2.92 ± 0.005 m and is 2.54 cm 
thick. It is mounted on a pair of spherical roller bearings 
held by a frame. The vessel top lid is installed in a 1.5 m 
diameter cylindrical flanged opening, and the inside lid 
surface is curved to complete the outer spherical bound- 
ary. The lid has four 13 cm diameter instrumentation 
ports centered at 60 cm cylindrical radius. Due to design 
constraints aimed at safe operation with liquid sodium 
metal as the working fluid, these four ports are the only 
penetrations through the outer boundary, and so are the 
only location from which we may make direct flow mea- 
surements. Port inserts hold measurement probes nearly 
flush with the inner surface of the outer sphere. 

The inner sphere has a diameter of 1.02±0.005 m and is 
supported on a 16.8 cm diameter shaft held coaxial with 
the outer shell by bearings at the bottom of the outer 
sphere and in the top lid. The inner sphere is driven from 
a 250 kW electrical motor through a calibrated Futek 
model TFF600 torsional load cell. The measured torque 
includes the torque from a pair of lip seals which add 
some confounding error. The outer sphere is driven by a 
250 kW induction motor mounted to the support frame. 
A timing belt reduction drive with a 25:3 ratio couples 
the outer sphere motor to a toothed pulley on the lid rim. 

Motor speeds are controlled to better than 0.2% by 
variable frequency drives, and optical encoders monitor 
the inner and outer sphere speeds. The drives estimate 



the motor torque from electrical current measurements 
and the torque estimate supplied by the inner motor drive 
agrees well with the calibrated strain-gauge torque sen- 
sor at motor speeds above about 2 Hz. This supports 
the use of the outer motor drive's reported torque as a 
reasonable estimate of the total torque required to drive 
the outer sphere, provided that the outer sphere angu- 
lar speed is above about 0.24 Hz, as it is for the data 
presented here. As the outer sphere speed is lowered be- 
low 0.24 Hz, the outer drive's reported torque becomes 
increasingly dominated by motor magnetizing current. 
The torque exerted by aerodynamic and bearing drag 
on the outer sphere is typically larger than the working 
fluid's contribution. However, this drag is repeatable and 
can be subtracted off by measuring the torque demanded 
with no differential rotation. 

A rotating computer acquires data from sensors in 
the instrumentation ports at a sampling rate of 512 Hz, 
recording data on a lab frame computer using a wire- 
less ethernet connection. Sensors include a Dantec model 
55R46 flush mount shear stress sensor driven by a TSI 
model 1750 constant temperature anemometer and three 
Kistler model 211B5 pressure transducers. The three 
pressure transducers are installed in three ports 90° apart 
on the 60 cm radius port circle. A thermocouple is used 
to monitor fluid temperature. 

A Met-Flow UVP-DUO pulsed Doppler ultrasound ve- 
locimeter is mounted in the rotating frame and paired 
with Signal Processing transducers, also communicating 
with the lab frame using wireless ethernet. Some velocity 
data in this paper was acquired with a Met-Flow UVP- 
Xl-PS ultrasound unit in the lab frame using a resonant 
transformer arrangement to couple the signal into the ro- 
tating frame. To scatter ultrasound, the flow is seeded 
with 150 /Ltm — 250 /zm polystyrene particles with nomi- 
nal density of 1.05 g/cm^. Limitations on the product of 
maximum measurable velocity and measurement depth 
constrain the velocity measurements in this paper to be 
very close to the wall, typically no more than 10 cm from 
the transducer face, which is either flush mounted with 
the wall or intruding no more than 10 cm into the flow. 
Intrusive transducers do not seem to measurably change 
the torque dynamics. When possible, we measure up- 
stream to avoid directly measuring the turbulent wake 
shed from an intrusive mounting scheme. 



III. TORQUE MEASUREMENTS 

The torques required to maintain constant speeds of 
the boundaries of in a rotating flow provide a global pic- 
ture of angular momentum transport and power dissipa- 
tion. Experimental torque measurements on wall-driven 
flows have been largely confined to the flow between con- 
centric cylinders^^^^^. The torque required to drive the 
flow between two spheres has received less attention, with 
previous results in the hydrodynamic turbulent regime 
seemingly limited to outer-sphere stationary measure- 



ments of liquid sodium^^^. Following Lathrop et al.^^, we 



1x10 



G 



1x10 




14x10 



10 
10 



Re 



10' 



10x10 



10 



G 



6x10 



2x10 



10 



(b) 



H ,. ^^ 



A 
▲ D 

LL 



-• -^ L 



*^0 0.5 1 1.5 2 2.5 3 3.5 

Ro 

FIG. 2. (a) The dimensionless torque G vs. the Reynolds 
number Re for the case flo ~ 0, stationary outer sphere, with 
a fit to G = ARe + BRe^ with A = 4.49 x lO" and B = 0.05. 
(fe) Mean value of the dimensionless torque G vs. Ro a.t E = 
2.1 X 10~^. Ro is varied by increasing the inner sphere speed 
in steps of Ro = 0.067, waiting 450 rotations per step. Three 
curves are differentiated by circles, triangles, and squares. In 
the ranges of Ro where the symbols overlap, flow exhibits 
bistable behavior. H and L denote the torque curves of the 
"high torque" and "low torque" states. There is a second 
bistable regime starting around Ro = 2.75, with LL labeling 
the lower torque state. 

define the dimensionless torque G on the inner sphere. 



G = 



T 



pv'^r-i 



(6) 



where p is the fluid density, v is the kinematic viscosity, 
ri is the inner sphere radius and T is the dimensional 
torque. The torque G as a function of Reynolds number 
Re, here deflned as 



Re = 



{Qi - flo)^^ 



(7) 



is bounded above by G oc Re^ from dimensional argu- 
ments. We flt the measured torque with the outer sphere 



stationary to G ex Re^ with a correction for linearly in- 
creasing seal drag. Data and fit with the outer sphere 
stationary are shown shown in Fig. 2(a) for the range of 
Re accessible in this experiment. At the lowest rotation 
rates, the measurement is significantly confounded by the 
bearing and seal torque. 

Figure 2(a) shows that G monotonically increases as 
Re is increased with the outer sphere stationary. This is 
not necessarily the case when we impart global rotation, 
provided that we do not hold Ro constant. When we 
hold the outer sphere speed constant (constant Ekman 
number) and super-rotate the inner sphere (Ro > 0) as 
in Fig. 2(6), we find that the mean torque on the inner 
sphere is non-monotonic as the dimensionless differential 
speed, Ro, and the Reynolds number. Re = Ro/ E are 
concurrently increased. As Ro is increased, the flow un- 
dergoes several transitions to significantly different tur- 
bulent flow states. In ranges of Ro in Fig. 2(6) where 
the torque data are presented in overlapping branches 
with different symbols, the flow exhibits bistable behav- 
ior with intermittent transitions between adjacent flow 
states. In these ranges of Ro, we plot the mean torque 
of each state, conditioned on state. 
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it in detail in this paper. 

A representative time series of the torque in the first 
bistable regime is shown in Fig. 3, with E ~ 2.1 x 10^^ 
and Ro — 2.13. The onset of the high torque state is 
abrupt, taking on the order of 10 rotations of the outer 
sphere. The torque overshoots the high state mean value 
at high torque onset by 10-15%. 

The end of the high torque state exhibits a slow decay 
of the torque to the low torque value, approximately ex- 
ponential with a time constant of 40 rotations of the outer 
sphere. In addition to the full transitions between the two 
torque levels, there are "excursions" where the torque de- 
cays toward the low mean value or rises toward the high 
mean value without fully reaching the other state. 

The qualitative bi-stability in Fig. 3 can be expressed 
quantitatively in the bimodal probability distribution of 
the torque shown in Fig. 4. The division of data into high 
torque and low torque states was done manually so as to 
exclude the transition regions. The resulting individual 
distributions of G for the high and low states based on 
this conditioning are shown in Fig. 4, as well as the full 
distribution. There is a small region of overlap in the 
H and L state individual distributions due to the differ- 
ence drawn between "transitions" and "excursions." The 
same manual division in states is used throughout the 
paper to condition other data on state. In Fig. 4, the 
high state mean torque is 1.4 times that in the low state. 
The torque fluctuations in the high torque state are con- 
siderably higher than in the low; the standard deviation 
of the high state torque is 1.8 times that in the low state. 
The low frequency fluid fluctuations responsible for this 
are also observed in the velocity and wall shear. 

In both Fig. 3 and Fig. 4 the torque data has been nu- 
merically low pass filtered. In Fig. 3 the cutoff frequency 
is 0.05 Hz. In Fig. 4, the filter cutoff frequency was cho- 
sen to be 0.5 Hz, where the fluid torque power spectrum 
appears to cross the mechanical vibration noise floor. In 
this way we retain the fastest measurable hydrodynami- 
cally relevant fluctuations. 



FIG. 3. Time series of G at fixed Ro = 2.13 und E = 
2.1 X 10~^, with time made dimensionless by the outer sphere 
rotation period. The raw torque signal has been numerically 
low pass filtered (/c — 0.05 Hz, 15 rotations of the outer 
sphere.). 

The state switching means that the transitions are 
not hysteretic, and long time mean torque through the 
switching range decreases with increasing Ro. We choose 
to plot the conditioned torque branches here instead of 
the long-time mean to emphasize our observation that 
these states seem the same as those present alone at 
higher and lower Ro. We will focus on the first bistable 
regime in this paper, between Ro — 1.8 and Ro = 2.75, 
and not the second bistable regime with lower state " LL" 
that begins near Ro — 2.75. This second bistable regime 
has several similarities to the first, but we will not discuss 



TABLE I. Statistics of the interval between high torque on- 
sets. At'ii is the time interval between two subsequent high 
torque onsets made dimensionless by ilo/2-K, so the time in- 
terval is measured in outer sphere rotations. Ro = 2.13, 
£; = 2.1 X IQ-^ 



(A4) 



C^At' 



Max(AfH) 



Min(AfH) 



717 



313 



1917 



390 



The interval between transitions is somewhat irregu- 
lar. Over 45 transitions at the parameters in Fig. 3, we 
observe the statistics shown in Table I. We also observe 
that the probability that the system is in one state or the 
other depends on Ro. Above a threshold value of Ro, we 
begin to observe state transitions to the low state, and 
the high torque state becomes less likely as Ro increases. 
Fig. 5 shows the probability that the system is in the 
high or low torque state for 4000 rotations of the outer 
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FIG. 4. Probability density of the dimensionless torque at 
Ro = 2.13, £ = 2.1 X 10"''. The full, unconditioned dis- 
tribution is denoted by small points. Solid circles denote 
conditioning on low torque state, and open circles on the 
high, with Gaussian solid and dashed curves for the low 
and the high respectively. The mean and standard devi- 
ation of the low torque state data are (G) — 6.82 x 10^" 
and ac ~ 2.74 x 10^. In the high torque state they are 
(G) = 9.53 x 10^" and oq = 5.02 x 10^ The data were low 
pass filtered at /c = 0.5 Hz to remove high frequency noise 
caused by mechanical vibration. 



sphere across the first bistable range of Ro. At values 
of Ro where transitions were not observed for more than 
4000 rotations of the outer sphere, a probability of one 
or zero was assigned. The probabihty that the system 
was in the low state was fit to 



Pr(L) 





1 — exp(— 7( 



Roo 



)) 



Ro < Roc 
Ro > Roc 



(8) 



with 7 = 8.25 and Roc — 1.80. The high torque prob- 
ability is given by Pr(iJ) = l-Pr(L). The physical im- 
phcation inherent in the exponential form is that there 
is no upper threshold where the high torque state be- 
comes impossible. Instead, it only becomes less likely as 
Ro is increased. However, the lower threshold for state 
transitions is well defined at Ron- 



IV. FLUID ANGULAR MOMENTUM 

The torque measurements presented so far only con- 
sidered the torque on the inner sphere. To see the com- 
plete picture of the angular momentum transport in the 
system, we examine the torque on both boundaries. The 
torque on the outer sphere is reported by the motor drive, 
though with less precision than the inner torque sensor 
measurement. 

The fluid filling the gap cannot undergo angular ac- 
celeration or deceleration for arbitrarily long times with 
the boundaries rotating at constant speeds. Therefore, 



FIG. 5. Probability that the flow was in the high torque (open 
circles) or low torque (closed circles) state as a function of 
Ro over the first bistable range with fixed iJ = 2.1 x 10"'. 
The dashed and solid lines are fits to the exponential form of 
Eq. (8). 



the time averaged net torque on the fluid must be zero, 
and the boundary torques must be equal and opposite, 
provided averaging is done over sufficiently long times. 
In the bistable regime, however, the net torque, shown in 
Fig. 6, reveals long periods of angular acceleration and 
deceleration interspersed with plateaus where the angu- 
lar momentum of the fluid remains nearly constant. The 
torques on the inner and outer boundaries only balance 
on averaging over many state transitions. We define the 
net torque on the fluid as 



T„, 



T + To 



(9) 



where T and To are the torques about the axis that the 
fluid exerts on the inner and outer spheres respectively. 
To is always negative when Ro > 0, as the fluid would 
tend to speed up the outer sphere rather than resist its 
motion. The total torque on the outer sphere as reported 
by the outer sphere motor drive, To^moas: is largely due 
to drag from the outer sphere bearings and air around 
the outer sphere. However, it is reproducible, and we 
deflne the torque that the interior fluid exerts on the 
outer sphere Tq in Eq. 9, as 



Tn = T 



o.mcas 



^dia 



(10) 



The drag, Tkrag is the outer drive torque required to keep 
the experiment in solid body rotation (with inner locked 
to outer) at the same outer sphere angular speed. In the 
Ro > flow states, some of the total torque required to 
maintain the speed of the outer shell against bearing and 
aerodynamic drag is supplied by the inner sphere via the 
fluid, and this is easily detected. 

The net torque is shown in Fig. 6 along with the sepa- 
rate inner and outer torques with the steady outer sphere 



bearing and aerodynamic drag subtracted. The dimen- 
sionless torques are defined as in Eq. 6, for example 
Gnct — T^nctl pv^Ti- Up to a constant of integration, we 
can calculate the fluid angular momentum L(i) about the 
rotation axis from the net torque, 



L{t) = \ T^etit')dt'. 



(11) 



Due to measurement limitations, we cannot integrate the 
torques from motor startup to fix the constant of inte- 
gration. Instead, we set the initial value of L to zero 
at an arbitrary time and make the resulting quantity di- 
mensionless by dividing it by the angular momentum the 
fluid would have if it were in solid body rotation with the 
outer sphere. 



L' = 



L 



-^fluidf^o 



(12) 



The moment of inertia of the fluid filling the gap is 
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(13) 



with a value of (1.14 ± 0.02) x lO"* kg m^. 

At the onset of the high torque state, as shown in 
Fig. 6(a), there is some prompt response of the outer 
sphere torque, indicating a certain amount of increased 
angular momentum transport. However, the increase in 
the torque on the outer sphere is insufficient to fully op- 
pose the increased inner sphere torque. At this point, 
the net torque becomes steadily positive, and the fluid 
accelerates. As this happens, the fluid torque Go on 
the outer sphere tends to slowly decrease in magnitude, 
though with large fluctuations. Eventually the inner 
sphere torque G starts the slower transition to the low 
torque state. At a point during the high to low transition, 
as shown in Fig. 6(6), the net torque becomes negative 
and the total angular momentum starts to decrease. The 
torque on the outer sphere continues a slow decay toward 
a value opposite and equal to that on the inner sphere, 
occasionally reaching a net torque fluctuating about zero 
as in times after Fig. 6(c). 

When both spheres rotate the angular momentum fluc- 
tuations in the range of Ro considered here may indi- 
cate a slow "store and release" process where long last- 
ing imbalances in the torques on the inner and outer 
spheres lead to fluid spin up and spin down. We note 
that the average magnitude of torque is similar between 
the outer stationary and the outer rotating cases for the 
same Reynolds number in Fig. 2(a) and Fig. 2(b). In 
both cases the mean torque on the inner sphere is be- 
tween G = 10^° and G = 10^^. Although the dynamics 
of the angular momentum transport seem quite different, 
the magnitude of the transport has not changed drasti- 
cally. 



FIG. 6. The angular momentum, inner torque, outer torque, 
and net torque. Ro — 2.13, E = 2.1 x 10~^. The upper plot 
shows the dimensionless angular momentum L' as defined in 
Eq. (12). The lower plot shows the inner torque G, the outer 
torque Go, and their sum Gnct. The bearing and aerodynamic 
drag on the outer sphere have been subtracted off, and the 
torques have been low pass filtered as in Fig. 3. 



V. MEAN FLOW MEASUREMENTS 

In the bistable regime, all measured flow quantities un- 
dergo transition along with the inner torque. Fig. 7 shows 
that the time averaged dimensionless wall shear r^ and 
measured dimensionless azimuthal velocity u' decrease 
sharply the onset of the high torque state. 

In Fig. 7, the velocity was measured in a shallow range 
near the surface of the outer sphere at 23.5° colatitude 
(60 cm cylindrical radius). This location is about 9 cm 
outside a vertical cylinder tangent to the inner sphere 
equator (called simply the "tangent cylinder" from here 
on). In Fig. 7, the transducer beam was in a plane normal 
to the cylindrical radius at the port and was inclined 
23.5° from pointing straight down (23.5° inclined from 
negative z). Since the transducer responds only to the 
velocity component along the beam axis, this orientation 
made it sensitive to the cylindrical radial, azimuthal, and 
vertical velocity components, Ug, u^, and u^. 

However, by using a remotely controlled rotatable 
mount, we determined that the time- averaged measured 
velocity is dominated by the azimuthal component, u^. 
When the transducer is rotated 90° from the usual orien- 
tation, making it least sensitive to u^, the mean velocity 
was about 10% of that seen when it maximally responds 
to the azimuthal flow. This indicates a meridional circu- 



lation of some importance. But it also m.eans that the 
meridional circulation provides a small contribution to 
the mean measured velocity when the transducer is ori- 
ented to be most sensitive to azimuthal flow. This is the 
case in Fig. 7, and so we treat the measured mean ve- 
locity as entirely azimuthal and correct for the vertical 
inclination, so that: 



{u^)t 



\Umca.s/t 

sin 23.5°' 



(14) 



We then divide by the outer sphere equatorial tangential 
velocity to make {u^)t dimcnsionless: 



fl r 



(15) 



This non-dimensionalization means that u' can be inter- 
preted as a local Rossby number. 

The wall shear stress sensor was calibrated against the 
measured torque T on the inner sphere with the outer 
sphere stationary. We assumed that the mean wall shear 
Tw on the outer sphere was 



T 



47rr3 cos 23.5°' 



(16) 



We fit the bridge voltage V and mean wall shear calcu- 
lated from the measured torque by Eq. (16) to 



V^ = At2/3 + BtI^' + C, 



(17) 



as was done by Lathrop et al.^^ We then used the cali- 
bration coefficients A,B and C to calculate the wall shear 
stress Tw from the measured bridge voltage for all subse- 
quent data with the outer sphere rotating. 

In the torque switching regimes, the slow fluctuations 
in wall shear and the mean azimuthal velocity shown in 
Fig. 7 are both strongly anti-correlated with the inner 
sphere torque. The simultaneous wall shear and torque 
data suggest that the two different states have different 
latitudinal distribution of the shear stress on the outer 
sphere. Fig. 6 shows that the torque the fluid exerts on 
the whole outer sphere is indeed somewhat greater in the 
high torque state. However, the mean shear stress at the 
measurement location in the high torque state is 65% of 
that measured in the low torque state (see also Fig. 8(6)), 
implying a shear stress concentration at high latitudes in 
the low torque state relative to the distribution of shear 
stress in the high torque state. 

These shear and velocity measurements suggest a fast 
central zonal flow in the low torque state that ceases sud- 
denly at the high torque onset. When the fluid around, 
above, and below the inner sphere is circulating faster 
in the low torque state, there is less drag on the in- 
ner sphere. The unusual aspect of this is that it takes 
less torque and less power input to maintain the faster 
circulation. This could indicate a transport barrier to 
energy and angular momentum in the low torque state. 
This change in transport is also important in interpret- 
ing the observation that the total angular momentum 
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FIG. 7. Simultaneous time series of velocity, wall shear stress, 
and torque. A space-time diagram of the low pass filtered 
velocity u (defined in Eq. 15) is shown at the top. This mea- 
surement is dominated by the azimuthal velocity u^. The 
velocity u is made dimensionless by the outer sphere tan- 
gential velocity, and so can be interpreted as a locally mea- 
sured Rossby number. The dimensionless wall shear stress 
t4 = 47rrorw/(pi^^r-i) is shown in the middle, and the di- 
mensionless torque G is shown at the bottom. The wall shear 
stress and torque have been low pass filtered with /c = 0.05iif z 
as before, which is comparable to the time averaging of the 
velocimetry. 



is often decreasing in the low torque state (See Fig. 6, 
line (6)), while the measured circulation is increasing or 
steady. This is only possible with a change in shape of 
the angular momentum profile as a function of cylindrical 
radius. An angular momentum transport barrier could 
explain this, possibly one associated with a fast zonal 
flow^''''^"''^'^^"^^. The anti-correlated torque and flow 
measurements will be discussed more in Sec. VIII. 



VI. TURBULENT FLOW FLUCTUATIONS 

The turbulent fluctuations are significantly different in 
the two flow states. The measured velocity fluctuations 
in the high torque state are larger than those in the low 
torque state, despite the lower mean velocity. Fig. 8(a) 
shows the probability density of the dimensionless veloc- 
ity conditioned on the low and high torque states, with 
Gaussian curves for comparison. The mean azimuthal 
velocity measured in the low torque state is 2.45 times 
that seen in the high torque state (in agreement with 
Fig. 7), while the fluctuations in the high torque state 
are 1.5 times the low torque state fluctuations. Based on 
this measurement, the high torque state has a consider- 
ably higher turbulence intensity (t„//(u'), of 57%. The 
low torque state turbulence intensity is 16%. 

The wall shear stress fluctuations are similar in magni- 
tude between the two states, as shown in Fig. 8(6). The 
low torque state mean wall shear is about 1.5 times that 
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FIG. 8. (a) The probability density of the dimensionless ve- 
locity conditioned on the torque state, with solid circles de- 
noting the low torque state and open circles denoting the high. 
Dashed and solid lines are Gaussian. Standard deviations are 
a"„' = 0.043 for the high state, cr„/ = 0.029 for the low state. 
The high state mean is (u') = 0.076, and the low state mean 
is (m') = 0.186. The means are dominated by azimuthal ve- 
locity, though the transducers are equally sensitive to u^ and 
M0 in this measurement. (6) The probability density of the 
dimensionless wall shear conditioned on torque. Solid circles 
again denote the low torque state, and open circles the high, 
with solid and dashed Gaussian curves. The mean and stan- 
dard deviation in the low torque state are (t^) — 1.63 x 10^^ 
and a^i — 1.81 x 10^". In the high torque state they are 
«) =1.06 X 10" and a^, = 1.71 x 10^". 



in the high torque state, and the standard deviations 
in the two states agree to within 5%. The high torque 
state wall shear distribution is significantly skewed, with 
a skewness of about 0.6. 



VII. WAVES 

Some of the most energetic fluctuations in the low 
torque state are system-scale wave motions with well de- 
fined frequencies. The high torque state has more broad- 
band, low frequency fluctuations. The conditional wall 



shear stress power spectra of Fig. 9 show this clearly. 
Two wave motions feature prominently in the low torque 
state but peaks are nearly absent in the high torque state. 
The lower peak in the low state spectrum of Fig. 9 has 
a frequency of 0.18 f2o- Measurements with the multi- 
ple pressure sensors show that this wave has azimuthal 
wavenumber m = \. The higher peak has a frequency of 
0.71 r^o at this Rossby number, and several higher har- 
monics are also visible. This higher frequency wave ap- 
pears to have to = 2. 

The frequency and azimuthal wavenumber of the lower 
frequency wave are consistent with the full sphere iner- 
tial mode'^ (3,1,-0.1766) in the {l,m,uj/^o) notation of 
Greenspan^ and Kelley et al7 . The inertial modes are 
not known analytically for the spherical shell'^"^, but 
there is good experimental agreement in frequencies and 
spatial patterns between inertial modes observed in tur- 
bulent spherical Couette flow^ for Ro < and modes 
of the full sphere^. Like the inertial modes observed 
previously^, the dependence of this mode's frequency on 
the Rossby number is weak. The spectrogram of Fig. 10 
shows the variation of pressure power spectra from a sin- 
gle pressure sensor as Ro is varied and E held constant. 
Line (a) in Fig. 10 is at 0.16 fio- The strong lower fre- 
quency peak starts higher, 0.18 f2o, and varies down to 
0.14 r^o and back up to 0.16 f^o as Ro is increased. This 
peak is featured in flow power spectra over a wide range, 
not disappearing until about Ro =12. It is worth men- 
tioning that (3, 1, —0.1766) inertial mode of the sphere is 
one of a special class of slow, geostrophic inertial modes 
that are equivalent to Rossby waves propagating on a 
solid body background^^ . 
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FIG. 9. Power spectra of wall shear stress at Ro = 2.13 and 
E — 2.1 x lO"'^, conditioned on state. Angular frequency 
has been made dimensionless using the outer sphere angular 
speed. The black curve is the spectrum from the low torque 
state, and the gray curve is that of the high torque state. The 
low torque spectrum has prominent peaks at uj/flo = 0.18, 
0.71, and harmonics. In the high torque state there are broad 
peaks at uu/Qo ~ 0.40 and 0.53. 
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FIG. 10. A spectrogram of wall pressure at 23.5° colatitude 
shows evidence of several flow transitions as Ro is varied. 
E = 1.2, X 10"'' and 0.1 <Ro< 4.4, waiting 430 rotations 
per step with steps of Ro — 0.1. Instead of averaging the 
power spectra over an entire step of Ro, there are 30 spectra 
per step in Ro, so some temporal evolution is visible at a 
given Ro. In the L state, there are two strong waves, the 
lower, at (a), varies only slightly in frequency with Ro. The 
higher frequency wave at (fe) varies more strongly with Ro, 
suggesting that advection by the mean flow is important in 
setting its frequency. 



The frequency of the stronger, higher frequency m = 2 
wave varies more strongly with the Rossby number. The 
equation for Fig. 10 line (b) is 



n. 



= 0.25i?o + 0.16. 



(18) 



This variation with Ro suggests a Rossby wave that is 
doppler shifted by advection. Rossby waves propagate 
where there is a gradient of potential vorticity, which 
in an isothermal, incompressible fluid is the quotient of 
fluid absolute vorticity and fluid column height. In the 
/?-plane approximation common for Rossby waves^'''"^^, 
the dimensionless frequency of the waves is given by 



— = fc— -^ 



(19) 



where k is the wavenumber and U is the velocity of a 
mean flow advecting the waves. The parameter /3 is re- 
lated to the background potential vorticity gradient on 
which the waves propagate. In our system, this gradient 
could result from the topographical effect of the slop- 
ing boundaries''''^^ or from a gradient of relative vorticity 
caused by the mean azimuthal flow profile. We observe 
that the zonal flow velocity U varies linearly with the 
Rossby number Ro in the low torque state. Outside the 
tangent cylinder with the inner sphere super-rotating, 
both the topographical contribution to /3 and the prob- 
able contribution from the mean relative vorticity would 



make it negative, consistent with the positive intercept of 
Eq. 18. Therefore Eq. 19 and Eq. 18 are indeed similar, 
and suggest an advected Rossby-type wave. 

At -E = 10~^ and Ro = 2.9, in the low torque state 
range, Guervilly and Cardin observe an tti = 2 Rossby- 
type wave"'^. They do not report an additional m — 1 
wave. However, if the physical mechanism that gives rise 
to the TO = 2 mode there is robust to lowering the Ek- 
man number from 10"'^ to 10"^, it could be that the 
m = 1 wave is excited by the m = 2 advected wave 
when the inertial mode damping is low enough. Reynolds 
stresses from nonlinear waves in rotating systems can 
transport angular momentum by driving strong zonal 
flows^'^"^'^'^^'™, and the strength of these zonal flows com- 
pared to the wave amplitude grows rapidly as the Ekman 
number is decreased^^. Therefore, the waves could play 
an important role in the angular momentum dynamics 
we observe. 



VIM. DISCUSSION 

One possibility for the state transitions discussed here 
is the formation and destruction of a fast zonal circula- 
tion at small cylindrical radius, with a resilient barrier to 
transport in the low torque state. Such barriers are com- 
mon features in rotating turbulent flows^^''^^"'*^'^^"^^'^^. 
The fast mean flow and strong wall shear stress at high 
latitude in the low torque state, along with the falling 
total angular momentum of Fig. 6 mean that the two 
states must have a different shape of the profile of an- 
gular momentum with cylindrical radius {s^fl{s)). The 
low torque state flow must favor, on average, fast circula- 
tion at the center of the experiment and slower at larger 
radius. A simple possibility would be a central zonal cir- 
culation bounded by a sharp shear layer modulated by 
the advected waves discussed in Sec. VII. This is shown 
schematically in Fig. 11(a). A collapse of this shear layer 
could cause the abrupt low state to high state transition, 
suddenly surrounding the inner sphere with low angu- 
lar velocity fluid and increasing the torque. The very 
short timescale of the L to H transition, 10 rotations of 
the outer sphere, suggests that some sort of rapid mix- 
ing is responsible. When the fluid is angular momentum 
is closer to well mixed, we expect a relaxation to a less 
steeply varying angular velocity profile, as in Fig. 11(b). 

In the low torque state schematic. Fig. 11(a), the az- 
imuthal flow outside the shear layer is slow, but not 
locked to the outer sphere. The angular momentum mea- 
surements of Fig. 6 suggest that the Ekman suction on 
the outer boundary drains angular momentum from the 
region outside the shear layer faster than the flux across 
the zonal flow boundary can replenish it until an equilib- 
rium is reached or a L to H transition occurs. 

At the transition from the high torque to low torque 
state, the slow decay of the inner torque and increase of 
the observed fluid velocity and wall shear suggest that 
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IX. DYNAMICAL BEHAVIOR 




(a)L 



(bJU 



FIG. 11. A sketch of two possible mean flow states. The 
low torque state is labeled L and the high labeled H. The 
low torque state at (a) is characterized by fast zonal circu- 
lation near the core of the experiment and large amplitude 
waves. In the high torque state, the velocity profile varies 
more gradually. The zonal circulation has been destroyed by 
mixing across the transport barrier, so the fluid near the inner 
sphere is slower and the torque higher. 



the transport barrier bounding the central zonal flow has 
re-formed, but the inner sphere must still provide an- 
gular momentum and energy to spin up the fast central 
flow. At this point, since the exterior fluid is now weakly 
coupled to the fast inner sphere, it starts to spin down. 

We mention finally that changes in mean profiles of 
angular momentum by vigorous mixing have been stud- 
ied in rotating thermal convection^^"^''^. Mixing by con- 
vectively driven fluctuations tends to flatten the angular 
momentum profile associated with solid body rotation, 
and this results in faster velocities in the rotating frame 
at high latitudes'^. Increasing vigor of convective mixing 
of the conserved angular momentum results in a switch 
from a prograde equatorial zonal jet to a retrograde one 
in the models of Aurnou et alJ'^. Advectively flattened 
angular momentum in the bulk of turbulent swirling flows 
is well established in Taylor-Couctte flow^^'^'^'^^ when the 
outer vessel is stationary, though the profile is not so flat 
in outer-stationary spherical Couette^^. 

Since the total angular momentum of the system is 
fluctuating, we must be careful not to use arguments 
based on strict angular momentum conservation. An- 
gular momentum is freely exchanged with the rotat- 
ing walls. An alternative explanation based on angular 
momentum mixing from vigorous turbulent fluctuations 
could involve a profile closer to the quadratic solid body 
profile in the high torque state and an increasingly flat- 
tened angular momentum profile in the low torque state, 
and would not involve the formation and destruction of 
a mixing barrier. It does, however, require more than 
a simple redistribution of angular momentum to explain 
the results of Fig. 6. 
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FIG. 12. Time delay embeddings of the slow velocity fluctua- 
tions at y' = 0.02 for three values of Ro below, in, and above 
the bistable range. Arrows in (6) show the direction taken 
by the transitions. The time delay Ai is 4.5 rotations of the 
outer sphere in all three cases. The dimensionless velocity as 
defined previously has been scaled by Ro, roughly collapsing 
the mean velocity and fluctuation levels of each state. These 
data have been low pass filtered with /c = 0.1 Hz, or a period 
of 7.5 rotations of the outer sphere. 

We also investigated the evolution in phase space of 
H-H-L transitions using low dimensional time delay em- 
bedding of low pass filtered velocity data. Figure 12 is 
a 2D embedding of a low pass filtered velocity time se- 
ries. The filter is a 4th order Butterworth with a cutoff 
frequency /c ~ 0.1 Hz. This corresponds to a period of 
7.5 rotations of the outer sphere. This velocity signal 
was plotted against the same signal 4.5 rotations prior. 
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Figure 12 (a) is Ro = 1.8, the lower threshold of the first 
bistable range. At (6), we have Ro = 2.13, the same as in 
Fig. 3, where the bistable switching is evident, and in (c), 
Ro = 2.4, above the range where spontaneous transitions 
are observed. 

The velocity is rescaled here by the outer sphere tan- 
gential speed to make it dimensionless and by Ro, the ex- 
pected dimensionless velocity scale relative to the outer 
sphere. Figure 12(5) could be interpreted as a hetero- 
clinic connection between two turbulent attractors, with 
connections between the high torque flow state in (a) and 
the low torque flow state in (c). The arrows in Fig 12(6) 
show the direction of flow in phase space as the system 
undergoes several transitions. Below the critical Rossby 
number in Eq. 8, we expect HoL connection to be bro- 
ken. Above the bistable range, we recall the probability 
of state as a function of Ro of Eq. 8 and Fig. 5 and expect 
that connections between attractors weaken but do not 
break. 



X. CONCLUSIONS 

The experimental results presented here reveal novel 
turbulent multiple stability and turbulent flow transi- 
tions in a previously unexplored parameter range of very 
high Reynolds number, rapidly rotating spherical Cou- 
ette flow. The transitions appear to involve the formation 
and destruction of zonal flow transport barriers and ex- 
hibit strong waves and unusual angular momentum trans- 
port. These results suggest spherical Couette as another 
straightforward laboratory testbed for studying turbu- 
lent multiple stability and could help lead to a better 
understanding of similar phenomena in natural systems. 

We are grateful to the NSF Geophysics program for 
funding and to the two anonymous reviewers for many 
helpful suggestions. 
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